Frozen light in periodic metamaterials 

Alex Figotin and Ilya Vitebskiy 



Abstract. Wave propagation in spatially periodic media, such as photonic 
crystals, can be qualitatively different from any uniform substance. The 
differences are particularly pronounced when the electromagnetic wavelength is 
comparable to the primitive translation of the periodic structure. In such a case, 
the periodic medium cannot be assigned any meaningful refractive index. Still, 
such features as negative refraction and/or opposite phase and group velocities 
for certain directions of light propagation can be found in almost any photonic 
crystal. The only reservation is that unlike hypothetical uniform left-handed 
media, photonic crystals are essentially anisotropic at frequency range of interest. 
Consider now a plane wave incident on a semi-infinite photonic crystal. One can 
assume, for instance, that in the case of positive refraction, the normal components 
of the group and the phase velocities of the transmitted Bloch wave have the same 
sign, while in the case of negative refraction, those components have opposite 
signs. What happens if the normal component of the transmitted wave group 
velocity vanishes? Let us call it a "zero-refraction" case. At first sight, zero 
normal component of the transmitted wave group velocity implies total refiection 
of the incident wave. But we demonstrate that total reflection is not the only 
possibility. Instead, the transmitted wave can appear in the form of an abnormal 
grazing mode with huge amplitude and nearly tangential group velocity. This 
spectacular phenomenon is extremely sensitive to the frequency and direction of 
propagation of the incident plane wave. These features can be very attractive 
in numerous applications, such as higher harmonic generation and wave mixing, 
light amplification and lasing, highly efficient superprizms, etc. 
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1. Introduction 

Wave propagation in spatially periodic media, such as photonic crystals, is 
qualitatively different from that of any uniform substance. The differences are 
particularly pronounced when the wavelength is comparable to the primitive 
translation L of the periodic structure P EJ |21 EI 111 El CI ■ The effects of strong spatial 
dispersion culminate when the group velocity of a traveling Bloch wave vanishes. One 
reason for this is that vanishing group velocity always implies a dramatic increase in 
density of modes at the respective frequency. In addition, vanishing group velocity 
also implies certain qualitative changes in the eigenmode structure, which can be 
accompanied by some spectacular effects in electromagnetic wave propagation. A 
particular example of the kind is the frozen mode regime associated with a dramatic 
enhancement of the wave transmitted to the periodic medium El ^3 ^1 ^] ■ In 
this paper we consider several different modifications of the frozen mode regime, each 
related to a specific singularity of the electromagnetic dispersion relation. For every 
case we present specific examples of periodic dielectric structures supporting such a 
phenomenon. For the most part, we restrict ourselves to periodic layered media, which 
are photonic crystals with one dimensional periodicity. 

We start with the simple case of a plane electromagnetic wave normally incident 
on a semi-infinite periodic layered structure. In the case of normal incidence, there are 
two different kinds of frozen mode regime. The first one is associated with a stationary 
inflection point of the dispersion relation in Fig. Qlb). In this case the incident 
plane wave of frequency is converted into the frozen mode with infinitesimal group 
velocity and very large diverging amplitude. A detailed description of this phenomenon 
can be found in our previous pubHcations on this subject (HI El dl ^1 IT^ . 
Qualitatively different kind of frozen mode regime occurs in the vicinity of a degenerate 
photonic band edge, shown in Fig. dc). In this latter case, the incident wave of 
frequency ujd is reflected back to space, but not before it creates a frozen mode with 
huge diverging amplitude. Not every periodic structure can support the frozen mode 
regime. Specifically, in the case of a normal incident, a unit cell of the periodic layered 
array must contain at least three layers, of which two must display a misaligned in- 
plane anisotropy, as shown in Fig. [21 In photonic crystals with three dimensional 
periodicity, the presence of anisotropic constitutive component is not necessary. 

In Section 3, we turn to the case of oblique wave propagation, which is of 
particular interest here. The frozen mode regime at oblique incidence can occur 
when the normal component of the transmitted wave group velocity vanishes, while 
its tangential component remains finite. In such a case, the transmitted wave is an 
abnormal grazing mode with a dramatically enhanced amplitude and nearly tangential 
energy flux. Remarkably, the photonic crystal reflectivity at this point can be very 
low, implying almost total conversion of the incident plane wave into the abnormal 
grazing mode. A significant advantage of this phenomenon is that it can occur in much 
simpler periodic structures, compared to those supporting the frozen mode regime at 
normal incidence. Examples are shown in Fig. |31andE] The presence of anisotropic 
layers is still required. 

Yet another interesting modification of the frozen mode regime is an abnormal 
subsurface wave in the vicinity of a degenerate photonic band edge. A regular 
surface wave usually decays exponentially with the distance from the surface in either 
direction. By contrast, the abnormal subsurface wave is extremely asymmetric. It 
does decay rapidly outside the photonic crystal. But inside the periodic medium, its 
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amplitude sharply increases, and reaches its maximum at a certain distance from the 
surface. Only after that the field amplitude begins a slow decay, as the distance from 
the surface further increases. The profile of such a subsurface wave is similar to that 
of the frozen mode in the vicinity of a degenerate band edge. This phenomenon is 
briefly discussed in section 4. 




wavenumber k wavenumber k wavenumber k 



Figure 1. Fragments of dispersion relations containing different stationary 
points: (a) a regular band edge (RBE) g, (b) a stationary inflection appoint 
(SfP) 0, (c) a degenerate band edge (DBF) d. 
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Figure 2. Periodic stack capable of supporting k — w diagram with degenerate 
band edge. A unit cell L includes three layers: two birefringent layers Ai and A2 
with misaligned in-plane anisotropy and one isotropic B layer. The misalignment 
angle i/j = ipi — ip2 between adjacent anisotropic layers Ai and A2 must be different 
from and tt/2. 



2. Frozen mode regime at normal incidence 

The essence of the frozen mode regime can be understood from the example of a plane 
monochromatic wave normally incident on a semi-infinite periodic layered structure, 
as shown in Fig. [S] An important requirement, though, is that some of the layers 
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Figure 3. Periodic layered structure with two layers A and S in a primitive cell 
L. The A layers (hatched) are anisotropic with one of the principle axes of the 
dielectric permittivity tensor making an oblique angle with the normal z to the 
layers {sxz 7^ 0). The B layers are isotropic. The x — z plane coincides with the 
mirror plane of the stack. 
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Figure 4. Periodic layered structure with two layers A and B in a unit cell L. 
The A layer has inplane anisotropy, while the B layer can be isotropic. 



must display a misaligned in-plane anisotropy. The simplest periodic layered array of 
the kind is shown in Fig. |21 

To start with, let us introduce some basic notations and definitions. Let ^'/, 
\E'fl, and be the incident, reflected and transmitted waves, respectively. Assume 
for now that all three monochromatic waves propagate along the z axis normal to 
the boundary of semi-infinite periodic layered structure in Figs. The photonic 
crystal boundary is located at z = 0. Electromagnetic field both inside (at z > 0) and 
outside (at z < 0) the periodic stack is independent of the x and y coordinates, and 
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its transverse components can be represented as a column- vector 



E^z) 

Ey{z) 
Hy{z) 



(1) 



where E {z) and H (z) are time-harmonic electric and magnetic fields. We also assume 
that anisotropic layers of the periodic array have in-plane anisotropy, so the fields E (z) 
and H [z) are normal to the direction of propagation 

E{z) L z,H{z) ^ z. (2) 

All four transverse field components in are continuous functions of z, which leads 
to the following standard boundary condition at z = 

^-T (0) = (0) -f (0) . (3) 

Note that the polarizations of the incident, reflected and transmitted waves can 
be different, because some of the layers of the periodic array display an in-plane 
anisotropy, as shown in the example in Fig. [3 The latter is essential for the possibility 
of the frozen mode regime. 

In periodic layered media, the electromagnetic eigenmodes 4'^ (z) are usually 
chosen in the Bloch form 

■^k {z + L) = e^'^'^^k (z) , (4) 

where the Bloch wavenumber k is defined up to a multiple of 2tt/L. The 
correspondence between u and k is referred to as the Bloch dispersion relation. Real 
k correspond to propagating (travehng) Bloch modes. Propagating modes belong to 
different spectral branches lo (k) separated by frequency gaps. The speed of a traveling 
wave in periodic medium is determined by the group velocity |3] 

u ~ dw/dk. (5) 

Normally, each spectral branch u (fc) develops stationary points ujg = uj (fc^) where the 
group velocity |SJ) of the corresponding propagating mode vanishes 

duj/dk = 0, at uj ^ LUs = uj (kg) . (6) 

Examples of different stationary points are shown in Fig. ^ where each of the 
frequencies LOg, too and ujd is associated with zero group velocity of the respective 
traveling wave. Stationary points © play essential role in the formation of frozen 
mode regime. 

By contrast, evanescent Bloch modes are characterized by complex wavenumbers 
A; = fc' -I- ik". Evanescent modes decay exponentially with the distance z from the 
boundary of semi-infinite periodic structure. Therefore, under normal circumstances, 
evanescent contribution to the transmitted wave ^'^ (z) can be significant only in close 
proximity of the surface. The situation can change dramatically when the frequency uj 
approaches one of the stationary point values ajg. At first sight, stationary points (|SJ| 
relate only to propagating Bloch modes. But in fact, in the vicinity of every stationary 
point frequency Wg , the imaginary part k" of the Bloch wavenumber of at least one of 
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the evanescent modes also vanishes. As a consequence, the respective evanescent mode 
decays very slowly, and its role may extend far beyond the photonic crystal boundary. 
In addition, in the special cases of interest, the electromagnetic field distribution^' {z) 
in the coexisting evanescent and propagating eigenmodes becomes very similar, as uj 
approaches lUs- This can result in spectacular resonance effects, such as the frozen 
mode regime. What exactly happens in the vicinity of a particular stationary point 
((nj essentially depends on its character and appears to be very different in each of the 
three cases presented in Fig. ^ In the next subsection we present a simple qualitative 
picture of the frozen mode regime based on energy conservation consideration. Then, 
in the following subsection, we provide more consistent analysis of the phenomenon. 



Semi-infinite 
photonic slab 




z = 



Figure 5. Plane wave normally incident on a semi-infinite photonic crystal. The 
subscripts /, R, and T refer to the incident, refiected and transmitted waves, 
respectively. In all cases, the amplitude of the incident wave is unity. 



2.1. Energy density and energy flux at frozen mode regime 

Let 5*7, Sr and St be the energy flux in the incident, reflected and transmitted waves 
in Fig. |S1 The transmission and reflection coefficients of a lossless semi-infinite medium 
are defined as 

St Sr ,„s 

T = TT, p = —5-' 

where 

Si + Sr = St, p = 1 - t. 

With certain reservations, the energy fiux St of the transmitted travelling wave can 
be expressed as 

St = Wtu, (8) 

where u is the group velocity, which is also the energy velocity, and Wt is the energy 
density 

Wt oc |*tP • 



Frozen light in periodic metamaterials 



7 



Evanescent modes do not contribute to the normal energy flux St in the case of 
a lossless semi-infinite periodic structure. Besides, evanescent contribution to the 
transmitted wave becomes negligible at a certain distance z from the surface. The 
simple expression © may not apply when the transmitted wave involves two or more 
propagating Bloch modes, but we will not deal with such a situation here. 

Vanishing group velocity u in JSJ implies that the transmitted wave energy flux 
St also vanishes, along with the respective transmission coefficient r in (jjj). The 
only exception could be if the energy density Wt of the transmitted wave increases 
dramatically in the vicinity of stationary point frequency. In other words, if Wt in 
((SJ grows fast enough, as lo approaches uJs, the product Wtu in ((SJ can remain finite 
even a,i lo — lo a- In such a case, a significant fraction of the incident radiation can 
be converted into the slow mode inside the semi-infinite periodic array. Hereinafter, 
the effect of a dramatic growth of the transmitted wave amplitude in the vicinity of 
a stationary point © will be referred to as the frozen mode regime. The possibility 
of such an effect is directly related to the character of a particular stationary point. 
From this point of view, let us consider three different situation presented in Fig. ^ 



2.1.1. Regular hand edge We start with the simplest case of a regular photonic band 
edge (RBE) shown in Fig. Q^a). It can be found in any periodic array, including any 
periodic layered structure. Just below the band edge frequency uog, the dispersion 
relation can be approximated by a quadratic parabola 

LO g — LJ <yi [k — kg f' , at w ^ Ug. (9) 

This yields the following frequency dependence of the propagating mode group velocity 
u below the photonic band edge 



— CX [kg — k) OC {Ujg — LU)^^^ , at CJ ^ UJg. (10) 

arv 



Due to the boundary condition ||2Jl, the amplitude of the transmitted propagating 
Bloch mode in this case remains finite and comparable to that of the incident wave. 
Therefore, the energy flux © associated with the transmitted slow mode vanishes, as 
Lv approaches ujg 

St = Wtu cc < ^ ^ \ ' (n) 
1^ 0, at Lj > ivg 

Formula Hll|) expresses the well-known fact that in the vicinity of a regular band edge, 
a lossless semi-infinite photonic crystal becomes totally reflective. 

2.1.2. Stationary inflection point A completely different situation occurs in the 
vicinity of a stationary inflection point (SIP) shown in Fig. H^b). Such a point can be 
found in periodic layered structures involving anisotropic and magnetic layers jHUHI, as 
well as in some photonic crystals with 2- and 3-dimensional periodicity. In the vicinity 
of a stationary inflection point wo? the dispersion relation can be approximated by a 
cubic parabola 

uj — bjQ (yi {k ~ kof" . (12) 

The propagating mode group velocity u vanishes as oj approaches loq from either 
direction ^ 

u — — (X (k — fco)^ (X (lu — ti^o)^^'^ ■ (13) 
ak 
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But remarkably, the amplitude of the transmitted propagating mode increases so that 
the respective energy density Wt diverges as uj ^ ujq 

M^T oc (tJ - ujoT"^'^ ■ (14) 

The expression (|13ll together with (|14|l yield that the energy flux of the transmitted 
slow mode remains finite even dX lo — ujq 

St = Wtu ^ Si, aX lu ujq. (15) 

The latter implies that the incident light is converted to the frozen mode with 
infinitesimal group velocity H13|l and huge diverging amplitude H14(l . Of course, in 
reality, the amplitude of the frozen mode will be limited by such factors as absorption, 
nonlinear effects, imperfection of the periodic dielectric array, deviation of the incident 
radiation from a perfect plane monochromatic wave, finiteness of the photonic crystal 
dimensions, etc. Still, with all these limitations in place, the frozen mode regime can 
be very attractive for a variety of practical applications. 

2.1.3. Degenerate band edge Let us turn to the case of a degenerate band edge (DBE) 
shown in Fig. ^c). Such a point can be found in dispersion relation of periodic layered 
structures involving misaligned anisotropic layers. An example is shown in Fig. [3 
Just below the degenerate band edge frequency cud, the dispersion relation lu (fc) can 
be approximated as 

ojd — uj oc {k — kd)'^ , ai Lo ^LOd (16) 
This yields the following frequency dependence of the propagating mode group velocity 

u ^ {kd - kf (X {ujd - t^)^^"* , at w ^ cjd. (17) 

dk 

The amplitude of the transmitted slow mode in this case diverges, as the frequency 
approaches the band edge value 

Wroclwrf-wj , aiuj-^ujDBE, (18) 

which constitutes the frozen mode regime. But the amplitude does not grow fast 
enough to offset the vanishing group velocity ((T7|l . Indeed, the expressions ((TTIl and 
H18|l together with ||SJ) yield for the energy flux 

a Tir f ('-^d — i^)^^^ , at to ^ LUd 

St = Wtu «i^ „\' ~ , (19) 
V, at uj > LUd 

implying that, in spite of the diverging energy density H18|) . the energy flux of the 
transmitted slow wave vanishes, as lu approaches LUd- 

The situation at a degenerate band edge H16|l can be viewed as intermediate 
between the frozen mode regime at a stationary inflection point l|12|l . and the vicinity 
of a regular band edge l^. Indeed, on the one hand, the incident wave at lu = LUd is 
totally reflected back to space, as is the case at a regular band edge. On the other 
hand, the transmitted field amplitude inside the periodic medium becomes very large 
as LU LUd, which is similar to what occurs at a stationary inflection point. The very 
large amplitude (|18|l of the transmitted wave at lu LUd still can be very attractive 
for applications. 

The above consideration does not explain the nature of the frozen mode regime, 
nor does it address the problem of the eigenmode composition of the frozen mode. All 
these questions are the subject of the next subsection. 



Frozen light in periodic metamaterials 



9 



2. 2. Physical nature of the frozen mode regime 

2.2.1. Eigenmode composition of transmitted wave In a periodic layered structure, 
at any given frequency lo, there are four electromagnetic eigenmodes with different 
polarizations and wavenumbers. But in the setting of Fig. |S1 where the semi-infinite 
periodic layered array occupies the half-space z > 0, the transmitted wave is a 
superposition of only two of the four Bloch eigenmodes. Indeed, the propagating 
waves with negative group velocity do not contribute to {z), nor do the evanescent 
modes exponentially growing with the distance z. Generally, one can distinguish three 
different possibilities. 

(i) Both Bloch components of the transmitted wave are propagating modes 

^T{z) = ^prl{z)+'^pr2{z), Z>Q. (20) 

The transmitted wave in (j2()|l is composed of two Bloch eigenmodes with two 
different real wavenumbers fci and fc2 and two different group velocities ui > 
and U2 > 0. This constitutes the phenomenon of double refraction, provided that 
ui and U2 are different. The other two Bloch components of the same frequency 
have negative group velocities and cannot contribute to the transmitted wave ^t- 

(ii) Both Bloch components of '^t are evanescent 

*T (^) = *e„l (^) + (2) , ^ > 0. (21) 

The respective two values of k are complex with positive imaginary parts k" > 0. 
This is the case when the frequency to falls into photonic band gap at uj > LUg 
in Fig. n^a) or at w > io^ in Fig. ^c). The fact that k" > implies that the 
wave amplitude decays with the distance z from the surface of the semi-infinite 
photonic crystal. In the case H21|l . the incident wave is totally reflected back to 
space by the semi-infinite periodic structure. 

(iii) One of the Bloch components of the transmitted wave is a propagating mode 
with M > 0, while the other is an evanescent mode with k" > 

(z) = (z) + ^lev (z), z> 0. (22) 

For example, this is the case at a; wq in Fig. H^b), as well as at w < tOg in Fig. 
n^a) and at oj < uj^ in Fig. ^c). As the distance z from the surface increases, 
the evanescent contribution '^ev in decays as exp {—zk"), and the resulting 
transmitted wave '^t (z) turns into a single propagating Bloch mode '^pr. 

Propagating modes with u < 0, as well as evanescent modes with k" < 0, never 
contribute to the transmitted wave '^t inside the semi-infinite stack in Fig. |S1 This 
statement is based on the following two assumptions: 

- The transmitted wave ^I't and the refiected wave ^'^ are originated from the 
plane wave incident on the semi- infinite photonic slab from the left, as shown 
in Fig. [S| 

- The layered array in Fig. ^occupies the entire half-space and is perfectly periodic 
at z > 0. 



Frozen light in periodic metamaterials 



10 



If either of the above conditions is violated, the electromagnetic field inside the 
periodic stack can be a superposition of four Bloch eigenmodes with either sign of 
the group velocity u of propagating contributions, or either sign of k" of evanescent 
contributions. For example, this would be the case if the periodic layered array in Fig. 
El had some kind of structural defects or a finite thickness. 

The propagating modes with u > and evanescent modes with k" > are referred 
to as forward waves. Only forward modes contribute to the transmitted wave \1/t (z) 
in the case of a periodic semi-infinite stack. The propagating modes with m < and 
evanescent modes with k" < are referred to as backward waves. Since the backward 
Bloch waves are not excited in the setting of Fig. [Sj they play no role in further 
consideration. 

Note that the assumption that the transmitted wave (z) is a superposition 
of propagating and/or evanescent Bloch eigenmodes may not apply if the frequency 
u) exactly coincides with one of the stationary point frequencies For example, at 
frequency ujq of stationary inflection point, there are no evanescent eigenmodes at all, 
and the transmitted wave ^'t (z) is a (non-Bloch) Floquet eigenmode linearly growing 
with z 0. The term non-Bloch means that the respective field distribution does not 
comply with the relation (@J) . Similar situation occurs at frequency tOd of degenerate 
band edge. At the same time, at any general frequency, including the vicinity of any 
stationary point ©, the transmitted wave ^'t (z) is a superposition of two Bloch 
eigenmodes, each of which is either propagating, or evanescent. 

In all three cases ()2()l - I22II , the contribution of a particular Bloch eigenmode to 
the transmitted wave ^'t depends on the polarization ^/ of the incident wave. One 
can always choose some special incident wave polarization, such that only one of the 
two forward Bloch modes is excited and the transmitted wave \E't is a single Bloch 
eigenmode. In the next subsection we will see that there is no frozen mode regime in 
the case of a single mode excitation. This fact relates to the very nature of the frozen 
mode regime. 

Knowing the eigenmode composition of the transmitted wave we can give a semi- 
qualitative description of what happens when the frequency co of the incident wave 
approaches one of the stationary points © in Fig. A consistent mathematical 
analysis of the asymptotic behavior in the vicinity of different stationary points is 
rather complicated. The details can be found in fT^ and references therein. 

2.2.2. Regular photonic band edge We start with the simplest case of a regular 
photonic band edge. There are two different possibilities in this case, but none of 
them is associated with the frozen mode regime. The first one relates to the trivial 
case where none of the layers of the periodic structure displays an in-plane anisotropy 
or gyrotropy. As the result, all Bloch eigenmodes are doubly degenerate with respect 
to polarization. A detailed description of this case can be found in textbooks on optics. 
Slightly different scenario occurs if some of the layers are anisotropic or gyrotropic and, 
as a result, the polarization degeneracy is lifted. Just below the band edge frequency 
ojg in Fig. n^a), the transmitted field ^'t (z) is a superposition of one propagating 
and one evanescent Bloch modes. Due to the boundary condition (PJ, the amplitude 
of the transmitted wave at z = is comparable to that of the incident wave. In the 
case of a generic polarization of the incident light, the amplitudes of the propagating 
and evanescent Bloch components at z = are also comparable to each other and to 



Frozen light in periodic metamaterials 



11 



the amplitude of the incident light 

IV (0)1 - |*e.„ (0)1 - |*/| , at < wg. (23) 

As the distance z from the surface increases, the evanescent component {z) 
decays rapidly, while the amplitude of the propagating component remains constant. 
Eventually, at a certain distance from the slab surface, the transmitted wave ^'t {z) 
becomes very close to the propagating mode 

{z) « ^'pr [z], at z > i, uj < ujg. (24) 

The evanescent component ev of the transmitted wave does not display any 
singularity at the band edge frequency LOg . The propagating mode "^pr does develop a 
singularity associated with vanishing group velocity at w — *■ — 0, but its amplitude 
remains finite and comparable to that of the incident wave. At w > ujg, this 
propagating mode turns into another evanescent mode in H21|) . The bottom line is that 
none of the Bloch components of the transmitted wave develops a large amplitude in 
the vicinity of a regular photonic band edge. There is no frozen mode regime in this 
case. 

2.2.3. Stationary inflection point A completely different situation develops in the 
vicinity of a stationary inflection point (|12|l of the dispersion relation. At w w Wq, the 
transmitted wave ^'^ is a superposition (|22|l of one propagating and one evanescent 
Bloch component. In contrast to the case of a regular photonic band edge, in the 
vicinity of both Bloch contributions to develop strong singularity. Specifically, 
as the frequency uj approaches ujo, both contributions grow dramatically, while 
remaining nearly equal and opposite in sign at the slab boundary [H] 

-^IpriO) ^ ~^ev{0) Oc\lU ~U!o\^^^^ , a.s LU LUq. (25) 

Due to the destructive interference l|25(l . the resulting field 

^-T (0) = V (0) + ^ev (0) 

at the photonic crystal surface is small enough to satisfy the boundary condition 
as illustrated in Figs. |H|and|7| As the distance z from the slab boundary increases, 
the destructive interference becomes less effective - in part because the evanescent 
contribution decays exponentially 

(2) « ^-e. (0) exp (-zfc") (26) 

while the amplitude of the propagating contribution remains constant and very 
large. Eventually, the transmitted wave '^t {z) reaches its large saturation value 
corresponding to its propagating component ^pr, as illustrated in Fig. [T] 

Note that the imaginary part k" of the evanescent mode wavenumber in H26|l also 
vanishes in the vicinity of stationary inflection point 

k" oc jw — woP^'^ , as — > luq, (27) 

reducing the rate of decay of the evanescent contribution H26|l . As a consequence, 
the resulting amplitude \1/t (z) of the transmitted wave reaches its large saturation 
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value '^pr only at a certain distance Z from the surface. This characteristic distance 
increases as the frequency approaches its critical value 

Z ocl/fc" oc |w - tJoT^^^ (28) 

If the frequency of the incident wave is exactly equal to the frozen mode frequency 
Wo, the transmitted wave {z) does not reduce to the sum (|22() of propagating and 
evanescent contributions, because at u; = wq there is no evanescent solutions to the 
Maxwell equations. Instead, "^t {z) corresponds to a non-Bloch Floquet eigenmode 
diverging linearly with z ^ 



(2) - *r (0) oc z^o, atw=a;o- 



(29) 
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Figure 6. Frequency dependence of squared amplitude of the transmitted wave 
(a), as well as its propagating (b) and evanescent (c) components at the surface of 
semi-infinite photonic crystal at 2 = 0. The stationary inflection point frequency 
is a;o . Due to the destructive interference 1^51 of the propagating and evanescent 
contributions, the amplitude (a) of the resulting transmitted field remains small 
enough to satisfy the boundary conditions. The amplitude of the incident wave 
is unity. 



2.2.4- Degenerate band edge While the situation with a regular photonic band edge 
((nj appears trivial, the case of a degenerate band edge (|16|l proves to be quite 
different. Below the degenerate band edge frequency ujd, the transmitted field is a 
superposition of one propagating and one evanescent components, while above 
LOd, the transmitted wave is a combination (|21l) of two evanescent components. In this 
respect, the regular and degenerate band edges are similar to each other. A crucial 
difference, though, is that in the vicinity of a degenerate band edge, both Bloch 
contributions to the transmitted wave diverge as u approaches lo^, which constitutes 
the frozen mode regime. 

Let us start with the transmission band at < w^- As the frequency to approaches 
(jjd — 0, both Bloch contributions in H22() grow sharply, while remaining nearly equal 
and opposite in sign at the surface 

(0) « -*et, (0) oc - wp^/"^, asw^Wd-O. (30) 
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Figure 7. Destructive interference of tlie propagating and evanescent components 
of the transmitted wave inside semi-infinite photonic crystal at z > 0. The 
frequency is close but not equal to that of the stationary inflection point, (a) The 
squared amplitude \^t{z)\^ of the resulting transmitted field - its amplitude 
at 2 = is small enough to satisfy the boundary conditions; (b) the squared 
amplitude \^pr (-2)1^ of the propagating contribution; (c) the squared amplitude 
l^'eu (^)P of the evanescent contribution. The amplitude of the incident wave is 
unity. The distance z from the surface is expressed in units of L. 



The destructive interference (|30|l ensures that the boundary condition © can be 
satisfied, while both Bloch contributions to {z) diverge. As the distance z from 
the slab boundary increases, the evanescent component ^'gy (z) dies out 

W « ^e. (0) exp (-zA:") (31) 

while the propagating component {z) remains very large. Eventually, as the 
distance z further increases, the transmitted wave '^t {z) reaches its large saturation 
value corresponding to its propagating component ^'p^ (z), as illustrated in Fig. |H1 
Note that the imaginary part k" of the evanescent mode wavenumber also vanishes in 
the vicinity of degenerate band edge 

k" cx |a; — o^til^^^ , as (jj ^ u)d^ (32) 

reducing the rate of decay of the evanescent contribution lp?T|l . As a consequence, 
the resulting amplitude '^t {z) of the transmitted wave reaches its large saturation 
value '^pr only at a certain distance Z from the surface. This characteristic distance 
increases as the frequency approaches its critical value LOd 

Zocl/fc"oc|w-Wd|"^/^ (33) 

as illustrated in Fig. llOf a. b). If the frequency uj of the incident wave is exactly 
equal to ojd, the transmitted wave '^t {z) does not reduce to the sum of two Bloch 
contributions. Instead, it corresponds to a non-Bloch Floquet eigenmode linearly 
diverging with z 

(z) - *T (0) oc z^d, at w = LOd- (34) 

This situation is illustrated in Fig. EJc). 

The above behavior appears to be very similar to that of the frozen mode regime at 
a stationary inflection point, shown in Fig. |7| Yet, there is a crucial difference between 
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the two cases. Indeed, according to (|19|) . in the immediate proximity of a degenerate 
band edge, the Pointing vector St of the transmitted wave '^t is infinitesimal, in spite 
of the diverging wave amphtude (|30|) . In other words, although the energy density 
Wt (X I^'tI^ of the frozen mode diverges as uj ujd — 0, it does not grow fast enough 
to offset the vanishing group velocity H17|l . As a consequence, the photonic crystal 
becomes totally reflective a.t oj = LUd- Of course, the total reflectivity persists at 
Lu > u)d, where there is no propagating modes at all. By contrast, in the case (|29ll of a 
stationary inflection point, the respective Pointing vector St is finite and can be even 
close to that of the incident wave. The latter implies low reflectivity and nearly total 
conversion of the incident wave energy into the frozen mode. 

Interesting situation occurs when we approach the degenerate band edge 
frequency from the band gap. In such a case, the transmitted field ^'t {z) is a 
superposition H21(l of two evanescent components. As the frequency uj approaches 
ujd- Both evanescent contributions grow sharply, while remaining nearly equal and 
opposite in sign at the photonic crystal boundary 

^'eul (0) « -1'e«2 (0) OC |cJd -t^r^^\ aS W ^ Wrf + 0. (35) 

Again, the destructive interference (|35(l ensures that the boundary condition is 
satisfied, while both evanescent contributions to ^'t {z) diverge in accordance with 
H35|l . As the distance z from the slab boundary increases, the destructive interference 
of these two evanescent components is lifted and the resulting field amplitude increases 
sharply, as shown in Fig. |51 But eventually, as the distance z further increases, 
the transmitted wave {z) completely decays, because both Bloch contributions to 
^T {z) are evanescent. The latter constitutes the major difference between the frozen 
mode regime above and below the DBE frequency LOd- The rate of the amplitude 
decay, as well as the position of the maximum of the transmitted wave amplitude in 
Fig. inia), are determined by the characteristic distance Z in (|33|l . 




20 40 60 20 40 60 20 40 60 

Distance z Distance z Distance z 

Figure 8. Smoothed filed distribution in the transmitted wave (a), along with its 
propagating (b) and evanescent (c) components at frequency slightly below that 
of the degenerate band edge in Fig. 0^c). The amplitude of the incident wave is 
unity. Similar graphs related to the stationary inflection point in Fig. 0^b) are 
shown in FiglSl 

2.2.5. Physical reason for the growing wave amplitude The wave ^'t {z) transmitted 
to the semi-infinite periodic layered medium is a superposition of two forward Bloch 
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Figure 9. Smoothed filed distribution in the transmitted wave (a), along with its 
two evanescent components (b) and (c), at a band gap frequency just above the 
degenerate band edge value u)^ in Fig. 0^c). Compare these graphs with those in 
Fig. Is] where the frequency lies in the transmission band below uJd- 



modes '^i (z) and 5*2 {z) 

*T {z) = *i {z) + *2 {z) . (36) 

The two eigenmodes in ()36|l can be a propagating and an evanescent, as in the case 
H22|l . or they can be both evanescent, as in the case H21|) . The physical reason for 
the sharp increase in the transmitted wave ampUtude at the frozen mode regime is 
that the two Bloch eigenmodes contributing to {z) become nearly indistinguishable 
from each other, as the frequency approaches its critical value 

^-1 (z) « a*2 [z) , as w ^ Ws, (37) 

where a is a scalar, and is the frozen mode frequency (wq or LOd)- 

Let us show under what circumstances the property 1)37(1 can lead to the frozen 
mode regime. The sum ((36|l of two nearly parallel column vectors and ^2 must 
match the boundary conditions © with the incident and reflected waves. If the 
incident wave polarization is general, then the nearly parallel Bloch components ^1 
and 4'2 must be very large and nearly equal and opposite 

*i (0) « -^-2 (0) , l^-i (0)1 « 1*2(0)1 » |*/|, (38) 

in order to satisfy the boundary conditions Indeed, since the incident field 

polarization is general, we have no reason to expect that the column vector \1/ (0) 
at the surface is nearly parallel to ^'i (0) and ^'2 (0). But on the other hand, the 
boundary conditions say that 

*(0) = *i(0) + *2(0) (39) 

Obviously, the only situation where the sum H39|l of two nearly parallel vectors can be 
not nearly parallel to either of them is the one described in ((38|) . 

There is one exception, though. Recall that in the vicinity of the frozen mode 
frequency, the two Bloch components ^'i and '^2 of the transmitted wave are nearly 
parallel to each other (see formula (|S3)- For this reason, if the polarization of the 
incident wave is such that \1/ (0) in H39|) is nearly parallel to any of the Bloch eigenmodes 
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Figure 10. Smoothed profile of the frozen mode at six difli^erent frequencies in the 
vicinity of degenerate band edge: (a) uj = u)ii — 10~'^c/L, (b) uj = uj^ — 10~^ c/L, 
(c) UJ = LUa, (d) uj = uja + lO'^c/L, (e) = oj^ + lO'^c/L, (f ) = cjd + lO'^c/L. 
In the transmission band at oj < tj^j, the asymptotic field value is very large and 
diverges as oj — > oj^. By contrast, in the band gap at ui > lUii, the asymptotic field 
value is zero. In all cases, the amplitude of the incident wave is unity. 

(0) and ^'2 (0), it is also nearly parallel to both of them. In this, and only this 
case, the amplitude of the transmitted wave ^'t {z) will be comparable to that of the 
incident wave. There is no frozen mode regime for this narrow range of the incident 
wave polarization. A particular case of the above situation is a regime of a single mode 
excitation, where only one of the two Bloch components or '^2 in H36|l contributes 
to the transmitted wave. 

Finally, let us reiterate that in the limiting cases oi to = ujq or uj — Ud, the 
transmitted wave "^/t {z) corresponds to a non-Bloch Floquet eigenmode H29|l or H34|) . 
respectively. Either of them linearly diverges with z. Again, the only exception is when 
the incident wave has the unique polarization, at which the transmitted wave {z) 
is a propagating Bloch eigenmode with zero group velocity and a limited amplitude, 
comparable to that of the incident wave. 
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3. Oblique propagation. Abnormal grazing modes 

A phenomenon similar to the frozen mode regime can also occur in the case of 
oblique wave propagation, where the incident, reflected and transmitted waves are 
all propagate at an angle to the z axis, as shown in Fig. 1111 Consider the situation 
where the normal component Uz of the group velocity of the transmitted propagating 
wave vanishes, while the tangential component remains finite. 



This is exactly what happens in the vicinity of the well-known phenomenon of total 
internal reflection (Sj. Similar phenomenon can be found in any photonic crystal 
at frequency corresponding to the transmission band edge for a particular direction 
of incidence. Remarkably, the total reflection of the incident wave is not the only 
possible outcome. Another alternative is that the transmitted wave forms an abnormal 
grazing mode with dramatically enhanced amplitude and tangential group velocity. 
The photonic crystal reflectivity at the critical point 1)40(1 can be even low - quite the 
opposite to what happens at total internal reflection. Low reflectivity implies that in 
spite of the vanishing normal component Uz of the transmitted wave group velocity, 
almost all the energy of the incident wave is converted into an abnormal grazing mode 
with drastically enhanced amplitude. The profile of such a grazing mode, i.e., the 
field dependence on the distance z from the surface, appears to be very similar to 
that of the frozen mode at normal incidence shown in Figs. |H1 El and E| The 
only difference is that the tangential component of the group velocity at an oblique 
version of the frozen mode regime can be large and comparable to the speed of light in 
vacuum. A significant advantage of the oblique version of the frozen mode regime is 
that it can occur in much simpler periodic structures, compared to those supporting 
the frozen mode regime at normal incidence. Examples of periodic layered arrays 
supporting only the oblique version of the frozen mode regime are shown in Figs. 
and ^ These structures are too simple to support any kind of frozen mode regime at 
normal incidence ~ they have only two different layers per unit cell, of which only one 
layer is anisotropic. But at oblique incidence, these relatively simple periodic arrays 
can support the frozen mode regime. The presence of at least one anisotropic layer in 
a unit cell L is still required. 

Further in this section we outline some basic properties of the frozen mode regime 
at oblique incidence. 

3.1. Axial dispersion relation: basic definitions 

Consider a plane monochromatic wave obliquely incident on a periodic semi-infinite 
stack, as shown in Fig. ^2 Due to the boundary conditions Q, the incident, reflected 
and transmitted waves should be assigned the same pair of tangential components 
kx, ky of the respective wave vectors 



while the axial (normal) components kz are all different. For the incident and reflected 
waves we have simply 




(40) 




(41) 




(42) 
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Figure 11. Scattering problem for a plane wave obliquely incident on a semi- 
infinite periodic layered medium. The arrows schematically shows the Pointing 
vectors of the incident, reflected and transmitted waves. The energy flux Sj of 
the incident plane wave is unity. 
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Figure 12. The case of oblique incidence, where the transmitted wave is a grazing 
mode with very large amplitude and purely tangential energy flux. 



Let us turn to the transmitted wave. The transmitted wave is usuaUy a 
composition of two Bloch eigenmodes with the same fcx = {kx,ky) from (|41|l . but 
different fc^. For given k± and w, the value of kz is obtained by solving the time- 
harmonic Maxwell equations in the periodic medium. The so found correspondence 
between the wavenumber k^ and the frequency lo at fixed k± is referred to as the axial 
or normal dispersion relation. Real k^ correspond to propagating (traveling) Bloch 
modes, while complex k^ correspond to evanescent modes, decaying with the distance 
z from from the surface. Unlike kx and ky, the Bloch wavenumber k^ is defined up to 
a multiple oi 2tt/L. 

Similarly to the case of normal propagation, the expression (|40|l defines stationary 
points of the axial dispersion relation. The definition H4U|I is a generalization of ® to 
the case of oblique propagation. Different kinds of axial stationary points are defined 
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as follows. 



A regular band edge of axial dispersion relation 

A stationary inflection point of axial dispersion relation 



dui 



0, 



afcf 



^0 



A degenerate band edge of axial dispersion relation 



du! 



0, 



5^ 



0, 



0, 



^0 



(43) 



(44) 



(45) 



The above definitions are analogous to those in lO, (|12|) . and (|16|) . related to 
the case of normal propagation. We still can refer to the k — u) diagrams in Fig. ^ 
where the quantity k is now understood as the normal component of the Bloch 
wavenumber at fixed fc_L. 

In order to further establish a close analogy between the cases of normal and 
oblique propagation, let us introduce some standard procedure commonly used in 
electrodynamics of layered media. In the case of oblique propagation, the time- 
harmonic electric and magnetic fields E [x, y, z) and H {x, y, z) depend on all three 
Cartesian coordinates, both inside (at z > 0) and outside (at z < 0) the periodic 
medium. Given the boundary conditions 141|l . the field dependence on the transverse 
coordinates x and y can be accounted for by the following substitution in the time- 
harmonic Maxwell equations 



E{r} = e*('=-^+'=«2')£ (z) , i? (r) = e^(''-^+'=«?')7i (z) 



(46) 



The next step is a separation of the transverse components of £ (z) and Ti, (z) into a 
close system of four linear differential equations 



d,^{z)=i-M{z)^{z), 
c 



(47) 



where ^' (z) is a vector-column 



*(z) 



£y{z) 

£x {z) 

£y{z) 



(48) 



similar to that in An explicit expression for the matrix M (z) can be found, for 
example, in Ref. along with the detailed discussion of its analytical properties. 

In the case of oblique incidence, the 4x4 matrix M (z) in the time-harmonic 
Maxwell equations H47|l depends not only on physical parameters of the periodic 
structure, but also on the tangential components fc_L = (fc^,, ky) of the wave vector. In 
the particular case of normal propagation, the symmetry of the 4x4 matrix M (z) in 
H47|l increases so that it may become incompatible with the existence of a stationary 
inflection point ifT^ and/or a degenerate band edge l(TB|l. This is why the axial frozen 



Frozen light in periodic metamaterials 



20 



mode regime at oblique incidence can occur in some periodic structures which are 
too simple to support the normal frozen mode regime. Examples such structures are 
presented in Figs. 13 and ^ 

All basic features of axially frozen mode regime at oblique incidence are virtually 
the same as in the case of normal propagation. In particular, all the formulae of 
the previous section describing the structure and composition of the transmitted field 
\1/t {z) remain unchanged. This close similarity holds both for the frozen mode regime 
at a stationary inflection point 1)44(1 and at a degenerate band edge (|45ll . In either 
case, Figs. 171 through [TUl give an adequate idea of the frozen mode structure. Still, 
there is one essential difference. Namely, in the case of axially frozen mode we have 
to remember that the tangential component of the group velocity is not zero, even if 
the normal component ((40|l vanishes. This means that the axially frozen mode is in 
fact an abnormal grazing mode with purely tangential energy flux, greatly enhanced 
amplitude, and very unusual profile. In the next subsection we elaborate on this point. 

3.2. Axially frozen mode as an abnormal grazing wave 

3.2.1. Stationary inflection point of axial dispersion relation Let us start with the 
case (|12|) of stationary inflection point. In the vicinity of luq, the transmitted wave 
V&T {z) is a superposition (|22|l of propagating and evanescent Bloch eigenmodes. 
According to H25|l . both contributions to (z) have huge and nearly equal and 
opposite values near the photonic crystal boundary, so that (0) is small enough 
to satisfy the boundary condition Q. As the distance z from the slab boundary 
increases, the evanescent component 5'e„ (z) decays exponentially, while the amplitude 
of the propagating component ^Pp^ (z) remains constant and very large. As the result, 
the field amplitude reaches its huge saturation value at a certain distance (|28|) from 
the slab boundary, as illustrated in Fig. [7| 

Let us now introduce a unit vector n parallel to the direction of incident wave 
propagation 

rix = kxc/uj, Uy = kyc/uj, = y/ 1 - (n^ + 7i2), (49) 

and let uq and uq be the direction of incidence and the frequency corresponding to a 
stationary infiection point H44() of the axial dispersion relation. When the frequency 
Lo tends to its critical value ujq at fixed direction of incidence uq, the saturation value 
vE'pr of the transmitted wave amplitude diverges as |a; — wqI '^^^ , according to H25|) . 
This behavior is illustrated in Fig. |S1 Conversely, when the direction of incidence ft 
tends to its critical value fto at fixed frequency wq, the respective saturation value "^pr 

_^ _^ 1/3 

of the transmitted wave amplitude diverges as |n — no| 

Consider now the frozen mode regime at oblique incidence in terms of refraction 
of the incident plane wave at the photonic crystal boundary. In Fig. E| we show a 
wide monochromatic beam of frequency u incident on the surface of the semi-infinite 
photonic crystal. The refraction angle 9t is determined by the ratio of the normal 
and tangential components of the group velocity 

— — Qrj. — arctan — . (50) 
2 

The direction of incidence no || 5/ in Fig. ^|is chosen so that the condition H44|) of 
the frozen mode regime is satisfied a,t uj — loq. As frequency lo tends to ljq from either 
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direction, the normal component of the transmitted wave group velocity vanishes, 
while the tangential component remains finite 

Uz ^ \lo — t^ol^'''^ ^ 0, ?r_L ^ uo as — > oJo- (51) 
This relation together with (|50|) yield 

— — oc — wqI" ^0 as uj ^ LOf). (52) 

Evidently, in the vicinity of the frozen mode frequency, the transmitted (refracted) 
propagating mode can be viewed as a grazing mode, because its group velocity becomes 
parallel to the surface. 

The most important and unique feature of this abnormal grazing mode directly 
relates to the fact that the transmittance of the photonic crystal remains finite even 
at a; = Wo- Indeed, let A/ and At be the cross-section areas of the incident and 
transmitted (refracted) beams in Fig. E| respectively. Obliviously, 

At _ cos 6>T , , 

Ai - cose I ^^'^> 

Let us also introduce the quantities 

Ui = AiSi, Ut = AtSt, (54) 

where Sj and St are the energy fluxes of the incident and transmitted waves. Uj and 
Ut are the total power transmitted by the incident and refracted beams, respectively. 
The expressions H53(l and H54|l imply that 

Ut _ St cos 0t _ {St)^ _ .j-^-n 
Ui " Sicos9i ~ (Si)^ ^ ' 

which is nothing more than a manifestation of the energy conservation. Equation (|55|l 
together with (|52|l yield 

c e I 1-2/3 

St — tSi — oc a; — Wo ^ oo as u; ^ ujq. (56) 

COSt^T 

where we have taken into account that r cos is of the order of unity, as a; — > . 

The expressions (|52l - I56f) show that in the vicinity of the frozen mode regime, 
the transmitted beam is a grazing mode with very large intensity St and sharply 
reduce cross-section area At, compared to those of the incident beam. By contrast, 
in the vicinity of the photonic band edge the transmittance r of the semi-infinite slab 
vanishes along with the normal component of the energy flux St of the transmitted 
(refracted) wave. 

The above qualitative analysis is only valid on the scales larger than the 
characteristic length Z in H28|l . The value Z defines the distance from photonic crystal 
surface, where the evanescent contribution to the transmitted wave is still significant. 
The latter restriction implies that the width of both the incident and the refracted 
beams must be much larger than Z. Otherwise, the transmitted wave cannot be 
treated as a beam, and the expressions (|52I ~ I56() do not apply. 



Frozen light in periodic metamaterials 



22 



Semi-infinite 
periodic sXack 




Figure 13. Incident and transmitted (refracted) waves in the vicinity of theaxial 
frozen mode regime regime. The reflected wave is not shown, dj and 9t are the 
incidence and refraction angles, Sj and St are the energy density fluxes of the 
incident and transmitted waves. Both the energy density and the energy density 
flux in the transmitted wave are much larger than the respective values in the 
incident wave. However, the total power transmitted by the refracted wave is 
smaller by factor r, due to much smaller cross-section area of the nearly grazing 
transmitted wave. 

3.2.2. Degenerate hand edge of axial dispersion relation Below the degenerate band 
edge frequency tOd in the transmitted wave '^t [z) is a superposition of 

propagating and evanescent Bloch eigenmodes. According to (I30|l . both contributions 
to \1/t {z) have very large and nearly equal and opposite values near the photonic 
crystal boundary at z = 0, so that '^t (0) is small enough to satisfy the boundary 
condition Q. As the distance z from the slab boundary increases, the evanescent 
component '^ev {z) decays exponentially, while the amplitude of the propagating 
component '^pr {z) remains constant and very large. As the result, the field amplitude 
reaches its huge saturation value at a certain distance from the slab boundary, as 
illustrated in Fig. |7| At first sight, this scenario is similar to the case of oblique frozen 
mode regime at stationary inflection point of the axial dispersion relation. But there is 
an important difference. As the frequency lo tends to its critical value w^;, the photonic 
crystal reflectivity approaches unity, implying total reflection of the incident radiation. 
In this respect, the situation below an axial degenerate band edge is reminiscent of 
the classical phenomenon of total internal reflection. The difference, though, is that 
the transmitted wave proflle '^t {z) is the same as in the previously discussed case of 
stationary inflection point of the axial dispersion relation. It implies that just below a 
degenerate band edge frequency of axial dispersion relation, the transmitted wave can 
be classifled as an abnormal grazing mode with large diverging amplitude and large 
tangential component of the energy flux. The typical profile of such a grazing mode 
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aX u) ^LOci is shown in Figs. Eta) and llUf a - b). 

If the incident wave frequency Hes sHghtly above Ud (in the band gap), the 
transmitted wave {z) is a superposition H21|) of two evanescent Bloch eigenmodes. 
The energy flux in the abnormal grazing mode now is strictly parallel to the surface, 
regardless of small variations in frequency and/or direction of propagation of the 
incident wave. This is a direct consequence of the total reflectivity of photonic crystals 
at band gap frequencies. The grazing mode profile in this case is shown in Figs. Efa) 
and llOf d - f). As the frequency oj approaches its critical value Wd, the amplitude 
and the thickness of the grazing mode diverge in accordance with (|35|l and 133|) . 
respectively. 

3.2.3. Surface waves in the vicinity of a degenerate band edge of axial dispersion 
relation So far in this section we have tacitly assumed that 



The inequality H57|l implies that the z component H42|l of the wave vector of the 
incident and reflected waves is real. This is a natural assumption when considering 
the problem of a plane wave incident on a semi-infinite photonic crystal. 
Consider now the opposite case where 



In this situation, there is no plane propagating waves in vacuum matching the 
boundary conditions H41|l at the surface. Still, if the frequency uj lies above the band 
edge of the axial dispersion relation (in the band gap), a superposition (|21|l of two 
forward evanescent eigenmodes with k'l > can produce a surface wave exponentially 
decaying with the distance z from the boundary in either direction. As the frequency 
uj approaches the degenerate band edge value ujd for a given fc_L, the surface wave 
profile can become similar to that of an abnormal grazing mode at ^ lj^;, shown 
in Figs. EI a) and llUf d - f). Although formally, it is still a surface wave, its profile 
now is highly unusual for a surface wave. Namely, the field amplitude in the periodic 
medium sharply increases with the distance z from the surface, reaches its maximum 
at a distance Z defined in H33|l . and only after that it begins a slow decay. 

4. Conclusion 

In this paper we outlined several different manifestations of the frozen mode regime 
in photonic crystals. Although all our numerical examples relate to periodic layered 
structures, in fact, the frozen mode regime is more like a universal wave phenomenon. 
Indeed, we can talk about different kinds of wave excitations in low-loss periodic 
media. But as soon as the respective Bloch dispersion relation displays a singularity 
like a stationary inflection point (|44|l or a degenerate band edge ()45|) . we have every 
reason to expect a very similar behavior involving the frozen mode regime. In other 
words, the possibility of the frozen mode regime is determined by some fundamental 
spectral properties of the periodic structure. If a periodic array is relatively simple, for 
instance, a stratified medium with one dimensional periodicity, its frequency spectrum 
may prove to be too simple to support the proper spectral singularity. The more 
complex the periodic structure is, the more likely it will be capable of supporting such 



kl + kl < lJc^. 



(57) 



kl + kl> <J-(?. 



(58) 
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Figure 14. Smoothed profile of the frozen mode in periodic layered structures 
composed of different number N of unit cells L in Fig. |^ The frequency is equal 
to that of degenerate band edge. The initial rate of growth of the frozen mode 
amplitude is independent of and described by Iij4l . The limiting case (d) of 
the semi-infinite structure is also shown in Fig. llUl cV In all cases, the incident 
wave has the same polarization and unity amplitude. 



a phenomenon. For instance, in the case of layered arrays we need birefringent layers 
and, sometimes, at least three layers in a unit cell. 

Another important question is how robust the frozen mode regime is. For 
instance, what happens if we introduce a small absorption or structural imperfections? 
Of course, these factors suppress the frozen mode amplitude. But in this respect, the 
frozen mode regime is no different from any other coherent or resonance effects it 
periodic structures. This problem can be solved at any particular frequency range by 
appropriate choice of the constitutive materials. 

Another fundamental restriction relates to the size of the periodic structure. In 
this paper we assumed that the periodic array occupies the entire half-space z > 0. 
A good insight on what happens to the frozen mode in a finite periodic array is given 
by Fig. 1141 These graphs prove that the frozen mode regime in a finite periodic array 
can be as robust as that in an imaginary semi-infinite structure. The optimal number 
of layers depends on such factors as the absorption characteristics of the constitutive 
materials, the geometrical imperfections of the periodic array, the desired degree of 
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field enhancement in the frozen mode, etc. On the other hand, in finite (bounded) 
photonic crystals, some new resonance phenomena can arise, such as transmission 
band edge resonances |21 El Qj- These effects, though, occur at slightly different 
frequencies and are qualitatively different from the frozen mode regime. 
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